Localization of electromagnetic waves in two-dimensional random dielectric systems 
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We rigorously calculate the propagation and scattering of electromagnetic waves by rectangular 
and random arrays of dielectric cylinders in a uniform medium. For regular arrays, the band 
structures are computed and complete bandgaps are discovered. For random arrays, the phenomenon 
of wave localization is investigated and compared in two scenarios: (1) wave propagating through 
the array of cylinders; this is the scenario which has been commonly considered in the literature, 
and (2) wave transmitted from a source located inside the ensemble. We show that within complete 
band gaps, results from the two scenarios are similar. Outside the gaps, however, there is a distinct 
fundamental difference, that is, waves can be blocked from propagation by disorders in the first 
scenario, but such an inhibition may not lead to inhibition or wave localization in the second scenario. 
The study suggests that the traditional method may be ambiguous in discerning localization effects. 
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X 



Wave localization is a peculiar property of random me- 
dia that completely block wave propagation due to mul- 
tiple scattering, thus inducing a surprising phase tran- 
sition, for example, in optical or acoustic transparency 
or electrical conductivity. When localized, waves remain 
confined in space until dissipated. 

More than two decades have passed since the phe- 
nomenon of wave localization was explored for propa- 
gation of electromagnetic (EM) waves in random media. 
During this period, a great body of literature has been 
generated 1 . And the interest in the subject continues to 
grow even further in recent years^^iSiLSiS. Despite the 
efforts, however, some important problems still remain 
unsolved. 

The first issue is with the way in which the localiza- 
tion effect is investigated. To date, claims of localiza- 
tion have been based on observations of the exponen- 
tial decay of waves as they propagate through disordered 
media™. That is, in most previous experimental or theo- 
retical studies, the apparatus was set up in such a manner 
that waves were transmitted at one end of a scattering 
ensemble, then the scattered waves were recorded on the 
other end to measure the transmission through the sam- 
ple. The results were then compared with the theory to 
infer the localization effect. In this method, it is quite 
plausible that other effects such as reflection and deflec- 
tion due to the presence of boundaries may also attenuate 
waves, resulting in a similar decay in transmission and 
obscuring the data interpretation. Therefore it is desir- 
able to look for a unique feature which can differentiate 
localization from other effects; the inability to discrimi- 
nate the localization effect from other effects has caused 
significant debate in the literature^. 

Second, although it has been suggested a while ago 
that the regions of localized states coincide with the po- 
sitions of the gaps, the relation between localization in 
random media and bandgaps of the corresponding regu- 
lar systems is still inconclusiv o 4 i 10 . 

Third, it has been the prevailing view over the past 
twenty years that all EM waves are localized in two 
dimensions (2D) for any given amount of disorder, fol- 
lowing the scaling analysis of electronic systemsAi. Re- 



cently, there is an intensive debate on this view from new 
experimentsASii^, as reviewed Since localization in 

electronic and EM systems has the same physical origin, 
it is therefore imperative to re-look at the view that all 
EM waves are always localized in 2D random systems. 
This task may be difficult, due to the obvious limitation 
of the finite sample size for either numerical or experi- 
mental workers, but at least one may examine whether 
the phenomenon of localization has been explored in a 
proper way in the past. 

With this Letter, we wish to shed new light on these 
questions. Here, we present a rigorous study of EM wave 
scattering and propagation in media containing many di- 
electric cylinders. The approach is based upon the self- 
consistent theory of multiple scattering^ and has been 
used previously to study acoustic localization in liquid 
media™ and acoustic attenuation by rigid cylinders in 
ai»i&. Wave propagation is expressed by a set of coupled 
exact equations and is solved rigorously. We show that 
wave localization can be achieved in ranges of frequencies, 
coincident with yet wider than the complete bandgap. 
For the phenomenon of wave localization, we compare 
two scenarios by analogy with the acoustic casei£: (1) 
the traditional setup of probing localization both numer- 
ically and experimnentally, as stated in, e. g. Ref4*£, 
that is, wave propagating through the array of cylinders, 
and (2) wave transmitted from a source located inside 
the ensemble. We show that within complete band gaps, 
results from the two scenarios are similar, whereas there 
is a distinct qualitative difference outside the gap. More- 
over, when localized, not only are waves confined near the 
transmitting source but a unique collective phenomenon 
emerges, illustrated by a phase diagram in analogy to the 
acoustic system^. 

The system considered here is similar to what has 
been presented in4. Assume that N uniform dielectric 
cylinders of radius a are placed in parallel in a uniform 
medium, perpendicular to the x — y plane. The arrange- 
ment can be either random or regular. For brevity, we 
only consider the case of the E-polarization, i. e. the E- 
field parallel to the z-direction. The qualitative features 
for both E- and H-polarizations are similar. The scatter- 
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ing and propagation of EM waves can be solved by using 
the exact formulation of Twerskjii^. While the details 
can be found ini&, here we brief the main procedures. A 
unit pulsating line source transmitting monochromatic 
waves is placed at a certain position. The scattered wave 
from each cylinder is a response to the total incident 
wave, which is composed of the direct contribution from 
the source and the multiply scattered waves from each 
of the other cylinders. The response function of a sin- 
gle cylinder is readily obtained in the form of the partial 
waves by invoking the usual boundary conditions across 
the cylinder surface. The total wave (E) at any space 
point is the sum of the direct wave (Eq) from the trans- 
mitting source and the scattered wave from all the cylin- 
ders. The normalized field is defined as T = E/Eq; thus 
the trivial geometrical spreading effect is eliminated. 
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FIG. 1: Left panel: The band structures computed by the 
plane wave expansion method. Right panel: Here is shown the 
normalized transmission log 10 \T\ 2 versus frequency; the solid 
line refers to the result from the [10] direction propagation, 
and the dotted line to that from the [11] direction propagation 
lines. 



the array of cylinders, labeled hereafter as the 'Outside' 
situation that imitates the traditional experimental^ and 
theoretical setups^, and (2) wave transmitted from a 
source located inside the ensemble, labeled hereafter as 
the 'Inside' situation. Both cases are illustrated by Fig. [2] 
For the 'Outside' case, all cylinders are randomly or regu- 
larly placed within a rectangular area with length L and 
width W. The transmitter and receiver are located at 
some distance from the two opposite sides of the scat- 
tering area. For the 'Inside' situation, all cylinders are 
placed either completely randomly or regularly within a 
circle of radius L with the source located at the center 
and the receiver located outside the scattering cloud. 
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FIG. 2: (a) The 'Outside' case: Electromagnetic propagation 
through a cloud of dielectric cylinders, (b) The 'Inside' case: 
Electromagnetic transmission from a line source located inside 
the array of dielectric cylinders. 



In line with4, the following parameters are used in the 
computation. The ratio of the dielectric constant be- 
tween the cylinders and the hosting medium is 10; the 
dielectric constant of the medium is taken as one. The 
filling factor f3, the fraction of area occupied by the cylin- 
ders per unit area, is 0.28. The radius a of the cylinders 
is 0.38 cm. The lattice constant d of the corresponding 
square latti ce ar ray of the cylinders is thus about 1.28 
cm (d = a^J-K I ft). For convenience, we scale all lengths 
by the lattice constant d. The computation is continued 
until the convergence is reached. 

First, in Fig. ^ we show the band structure of the cor- 
responding square lattice arrangement of the cylinders, 
obtained by the plane wave method. The wave trans- 
mission in two symmetric directions is also shown. Two 
complete bandgap regions are identified and are consis- 
tent with the highly attenuated regions in the transmis- 
sion computation. These results are also consistent with 
that in Fig. 4 o^, thereby verifying our numerical scheme. 

To investigate wave localization, two situations are 
considered and compared: (1) wave propagating through 



The frequency response of the averaged logarithmic 
transmission is presented in Fig. [3] for both 'Inside' and 
'Outside' scenarios. Here we see that the disorder some- 
what tends to enhance transmission within the bandgaps 
for both scenarios, while obviously reduces the transmis- 
sion for all frequencies outside the gaps in the 'Outside' 
case. For the 'Inside' situation, however, the reduction 
for regions outside the gaps is not generally obvious, and 
is only seen near the gap edges. It has been suggested 
in the literature that the transmission reduction in the 
'Outside' scenario indicates wave localization. For ex- 
ample, the authors in4<& computed the transmission in 
the context of the 'Outside' scenario, and subsequently 
obtained the localization length for all frequencies. We 
find that this approach towards localization may not be 
appropriate. The reasons follow. 

If the transmission reduction in the 'Outside' scenario 
is only caused by the localization effect, it will be implied 
that the random system only supports localized states. 
Then waves will not be allowed to propagate not only 
through but also inside the system. Therefore we would 
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expect the transmission to follow an exponential decay 
with increasing sample size for both 'Inside' and 'Outside' 
setups. 
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FIG. 3: Transmission versus frequency for both random and 
regular arrays of cylinders: (a) the 'Outside' case with W = 6 
and L = 10; and (b) the 'Inside' case with L = 10. Please refer 
to Fig. |5|and the text for the explanation about the 'Outside' 
and 'Inside' cases. 



Fig.EIpresents the results for the random ensemble av- 
eraged transmission and its fluctuation as a function of 
the sample size at two frequencies. The sample size is 
varied by adjusting the number of the cylinders. For the 
'Outside' case, we have done the following to remove the 
effect of the width W. With a fixed sample size (i. e. the 
length L), we plot the transmission versus width. We 
find that the transmission is very nicely saturating to a 
certain value in an exponential manner. We have done 
for several lengths, and obtained the corresponding sat- 
urated value for each length. Then we plot these values 
versus sample lengths. As an example, the results for 
8.64 GHz are shown in Fig. EI (e) and (f). For 6.54 GHz, 
the localization is strong, the width effect diminishes very 
quickly when the width increases. Here the plot for 6.54 
GHz has width 26 in the 'Outside' scenario. Note that 
the width should not be started at a value too close to 
zero; otherwise the variance will be too large, making 
the results unstable. The average has been taken for 500 
configuration to ensure the stability. 

A few important features are discovered. For the fre- 
quency of 6.54 GHz (within the first gap), the transmis- 
sion decays exponentially with the sample size for both 
'Outside' and 'Inside' situations with almost the same 
slop of -1.35, suggesting that at this frequency, waves 
are localized. And inside the localization regime, the ab- 
solute value of the transmission fluctuation is small, as 
expected from an earlier worki Here we see that within 
the localization regime, wave localization can be indeed 
observed in both 'Outside' and 'Inside' scenarios. 

For 8.64 GHz (between the first and the second gaps), 
the 'Outside' and 'Inside' scenarios differ significantly. 
For the 'Outside' case the transmission decreases expo- 
nentially with a slop of -0.0612 along the path. If this ex- 
ponential decay is caused by localization, then we should 
also observe the exponential decay for the same sample 
size (L) in the 'Inside' scenario. The result in the cen- 
ter panel of Fig. 0] clearly does not support this point of 



view. Instead, Fig. El tends to indicate that waves are 
not yet localized at 8.64 GHz in the 'Inside' scenario. 
The fact that the exponential decay only occurs in one 
scenario but not in the other for the same sample size 
(L) is itself intriguing and important. Therefore we may 
conclude that the 'Outside' scenario is inappropriate in 
isolating the localization effect, and it would be a mis- 
take to interpret the exponential decay or transmission 
reduction shown in the 'Outside' situation as a conclu- 
sive indication of wave localization. Furthermore, as at 
this frequency, waves are not yet localized in the 'Inside' 
case and they have a weaker exponential decay in the 
'Outside' case, the transmission will be more sensitive 
to the arrangement of the cylinders. Therefore the fluc- 
tuation at this frequency is stronger than that at 6.54 
GHz. However, the ratio between the fluctuation and 
the transmission at 8.64 GHz can be smaller than that 
at 6.54 GHz. 
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FIG. 4: The averaged logarithmic transmission and its fluctu- 
ation versus the sample size for two frequencies: one is within 
the first bandgap and the other is above the first but below 
the second gap. The left and center panels refer to the 'Out- 
side' and 'Inside' cases respectively. The estimated slops for 
the transmission are indicated in the figure. The right panel 
shows the effect of width W and the plot of the transmission 
versus length L at the extrapolated infinite width (see the 
text) . 



To this end, a few notes are appropriate. We have 
also examined other frequencies in general and two in 
particular: one is within the second gap and the other 
is above the second gap. The results are very similar to 
that shown in Fig. ^ For brevity, we will not show the 
results here. From Fig. |21 the fact that the transmission 
reduction occurs not only within but also outside the gaps 
(at areas around the edges of the gaps) indicates that 
the localized regions are coincident with the complete 
bandgaps, and these regions seem wider than the gaps. 
Our results show that although the disorders may block 
waves from propagation through the medium, but they 
may not yet localize the waves inside a 2D system. 

Now we discuss a unique feature of EM wave localiza- 
tion. The energy flow of EM waves is J ~ E x H. By 
invoking the Maxwell equations to relate the electrical 
and magnetic fields, we can derive that the time aver- 
aged energy flow is < J > t = h dtJ ~ \E\ 2 V9, where 
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the electrical field is written as E = e.E\E\e %e , with &e 
denoting the direction, \E\ and 6 being the amplitude 
and the phase respectively. It is clear that when 9 is 
constant, at least by spatial domains, while \E\ ^ 0, the 
flow would come to a stop and the energy will be lo- 
calized or stored in the space. We assign a unit phase 
vector, u = cos 9ie x + sin^e^ to the oscillation phase 0i 
of the dipoles. Here e x and e v are unit vectors in the 
x and y directions respectively. These phase vectors are 
represented by a phase diagram in the x — y plane. 
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FIG. 5: The phase diagram and spatial distribution of elec- 
tromagnetic energy for two frequencies for one random config- 
uration. Left panel: the phase diagram for the phase vectors 
defined in the text; here the phase of the direct field Eo is set 
to zero. Right panel: the energy spatial distribution. 



In Fig. the two-dimensional spatial distribution of 
EM energy (~ \E/Eq\ 2 ) and the phase vectors of the 
E-field are plotted for the two frequencies discussed in 
Fig. The phase vectors are located randomly in the 



x — y plane but to avoid the positions of the cylinders. 
The 'Inside' scenario is considered. Here we clearly see 
that for 6.54 GHz, the energy is mainly confined near the 
source, consistent with Fig.0] The phase vectors are or- 
derly oriented. These fully comply with the above general 
discussion. Therefore at this frequency, EM wave is in- 
deed localized. When we increasingly add an imaginary 
part to the dielectric constant, the ordered orientation 
of the phase vectors will disappear, confirming that the 
phase coherence is a unique feature of EM wave localiza- 
tion. We note from Fig.[S]that near the sample boundary, 
the phase vectors start to point to different directions. 
This is because the numerical simulation is carried out 
for a finite sample size. For a finite system, the energy 
can leak out at the boundary, resulting in disappearance 
of the phase coherence. When enlarging the sample size, 
we observe that the area showing the perfect phase co- 
herence will increase. At 8.64 GHz, however, there is no 
ordering in the phase vectors u(r). The phase vectors 
point to various directions. The energy distribution is 
extended in the x — y plane, and no EM wave localiza- 
tion appears, in agreement with what has been described 
for Fig. H 

In summary, we have examined some fundamental 
problems of EM wave localization in 2D. Although it may 
be still hard to conclude that extended waves are possi- 
ble in 2D random media, as limited by the finite sample 
size, the present results do indicate that the traditional 
method may be unable to isolate the localization effect. 
It is also shown that the localization region is related 
to and seems to be wider than the complete bandgaps. 
When localized, not only are waves confined near the 
transmitting source but a unique collective phenomenon 
emerges. 
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